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Abstract 


We study the capture cross-section of charged particles by a weakly charged Schwarzschild black hole. The 
dependence of the maximum impact parameter for capture on the particle’s energy is investigated numerically for 
different values of the electromagnetic coupling strength between the particle and the black hole. The capture 
cross-section is then calculated. We show that the capture cross-section is independent of the electromagnetic 
coupling for ultra-relativistic particles. The astrophysical implications of our results are discussed. 
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1. Introduction 


Studying the capture cross-section of black holes is central to 
understanding how matter interacts with them. It helps us 
understand the process of matter accretion by a black hole 
which in turn determines how its mass, angular momentum and 
charge evolve. It can also help us understand the environment 
near black holes. Moreover, scrutinizing capture cross-section 
can be used to test theories of gravity in strong gravitational 
fields. 

Astrophysicists generally assume black holes are electrically 
neutral. This is because they would quickly attract oppositely 
charged matter to balance out any access charge. However, 
there are compelling reasons why weakly charged black holes 
might exist as discussed in Al Zahrani (2021, 2022), Zajaéek & 
Tursunov (1904), Zajaček et al. (2019, 2018), Carter (1973) 
and the references therein. The differences in how a black hole 
accretes electrons and protons within its plasma environment, 
influenced by radiation, could render it charged. Also, the spin 
of a black hole in the presence of a magnetic field can induce 
the accretion of charged particles. In fact, using the Event 
Horizon Telescope (EHT) observations, it was inferred that Sgr 
A* and M87 can be charged (Ghosh & Afrin 2023; 
Kocherlakota et al. 2021). The black hole’s charge is weak in 
the sense that it has no tangible effect on spacetime, but its 
effect on charged particle dynamics is prominent. 

There are numerous astrophysical scenarios wherein charged 
particles are drawn into black holes. Stars within the Roche 
limits near black holes often contribute matter through tidal 
interactions. Additionally, stars emit streams of charged 
particles as stellar winds. Highly energetic charged particles, 
resulting from supernovae, gamma-ray bursts, and bipolar jets 
from compact objects, frequently find their way into the 
vicinity of black holes. These processes collectively enrich the 


environment around black holes with a significant population 
of charged particles. 

The concept of capture cross-sections has been explored 
extensively for various black hole types. Foundational 
treatment which examines photon and neutral particle capture 
by Schwarzschild black holes was given in several mono- 
graphs, such as Frolov & Zelnikov (2011). Further work 
addressed capture cross-sections of charged and neutral 
particles by Kerr-Newman black holes, including the 
implications for black hole spin and charge evolution (Young 
1976). Capture by Reissner—Nordstr6m black holes was also 
investigated (Zakharov 1994). In the context of higher- 
dimensional black holes, studies have focused on calculating 
photon critical impact parameters for Schwarzschild—Tangher- 
lini black holes (Connell & Frolov 2008; Tsukamoto et al. 
2014; Singh & Ghosh 2018; Bugden 2020). The capture cross- 
section for massive particles was determined in Ahmedov et al. 
(2021). Additionally, research extends to particle capture in 
Myers—Perry rotating spacetime which describes rotating black 
holes in five-dimensions (Gooding & Frolov 2008). Moreover, 
wave capture cross-sections have been studied for various 
black hole configurations (see Anacleto et al. 2023 and the 
references within). 

In this research, we examine the capture cross-section of 
charged particles by a weakly charged Schwarzschild black 
hole and discuss the astrophysical consequences of our 
findings. The paper is organized as follows: In Section 2, we 
review the dynamics of charged particles in the background of 
a weakly charged black hole. We then review the capture cross- 
section of neutral particles in Section 3. The capture cross- 
section of charged particles is calculated for different coupling 
strengths and particle energies in Section 4. Finally, we 
summarize our main findings and discuss their astrophysical 
consequences in Section 5. We use the sign conventions 
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adopted in Misner et al. (1973) and geometrized units where 
c=G=k=1, where k is the electrostatic constant. 


2. Charged Particles Near a Weakly Charged 
Schwarzschild Black Hole 


Here, we review the dynamics of charged particles near a 
weakly charged black hole. The spacetime geometry around a 
black hole of mass M and charge Q is described by the 
Reissner—Nordstr6n metric which reads (Misner et al. 1973) 


ds? = —hdt? + h-'dr? + r2d@? + r? sin? @dg?, (1) 


where h = 1 — rs/r + Q? /r and rs = 2M is the Schwarzschild 
radius. The electromagnetic 4-potential is 


A = 250. (2) 


oor 
However, when the charge is weak we can ignore the curvature 
due to it and use the Schwarzschild metric, which reads 
(Misner et al. 1973) 
ds? = —fdt? + f-'dr? + rd? + r° sin? 6dd?, (3) 


where f= 1 — rg/r. This weak charge approximation is valid 
unless the charge creates curvature comparable to that due to 
the black hole’s mass. This happens when 


Q? ~ M?. (4) 
In conventional units, the weak charge approximation fails 
when 
G'?M 
k! /2 


AJ 102° M 


coulombs. (5) 
© 

This charge is way greater than the greatest estimated charge on 
any black hole. Although the black hole charge is tiny, its effect 
on charged particle dynamics is profound because it is 
multiplied by the charge-to-mass ratio of these particles 
(~107! m™! for electrons and ~10!8 m™! for protons). 

The Lagrangian describing a charged particle of charge q and 
mass m in a spacetime described by a metric g,,, and an 
electromagnetic field produced by a 4-potential A” reads 
(Chandrasekhar 1983) 


1 
a al a qu"A,, (6) 


where u” = dx"”/dr is the particle’s 4-velocity and 7 is its 
proper time. In our case, the Lagrangian becomes 


PORO 


— qQ—. (7) 
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This Lagrangian is cyclic in t and ¢, which means that the 
particle’s energy and azimuthal angular momentum are 
constants of motion. The specific energy and azimuthal angular 
momentum are, respectively, given by 


1 OL dt 
aa malt) "ar i 
dr 

f= ae = r? sin? ge. (9) 

m a(2) dr 
Combining these equations with the normalization condition 

8 U4"u" = —1 and solving for dr/dr give 

2 2 2 2 
OES EE 8 +1] 
dr mr dr r? sin? 0 


(10) 


In the equatorial plane where 8 = 7/2, the equation becomes 


-e-gi-n ow 


Let us rewrite the last equation in a dimensionless form. We 
first introduce the following dimensionless quantities: 


, kae Tea. (12) 


7 
rs rs rs 


Equation (11) then becomes 
2 2 2 
(2) =[e A E (13) 
dT p p’ 


ga 22. (14) 


where 


The parameter a represents the relative strength of the 
electromagnetic force to the Newtonian gravitational force. 
We can rewrite Equation (13) as 


d 2 
(2) sVe =W: (15) 
where 
2 
W= 2 + ite n ‘|, (16) 


is an effective potential. It is V} that corresponds to physical, 
future-directed motion and hence will be used in all of the 
analyses below. Without loss of generality, we will consider 
£ > 0 only. 

It was estimated in Zajaéek et al. (2019) that the charge of 
Sgr A* is 10°10'° coulomb. Using the lower limit of charge, 
the coupling constants for electrons a, and protons a, near Sgr 


L 
A*, which has a mass of M=4.3 x 10°Mo according to 
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GRAVITY Collaboration (2023), are 
Qe ~ 10°, (17) 
Qp ~ 106. (18) 


3. Capture Cross-section of Neutral Particles 


Before we tackle the main problem, let us find the capture 
cross-cross section for neutral particles first. Setting a = 0, the 
effective potential V, reduces to 


2 
Ves (5 $ i) (19) 
p 


Capture occurs whenever the particle’s energy is greater than 
the maximum of V}. The function V, is at an extremum when 
dV, /dp=0 or 


pP? + (3 — 2p)L? = 0, (20) 
which gives the position of the extrema in terms of £ as 


p = P+ til —3, (21) 


where £L € [J3, co). When L = V3 (Lyin), p+ and p_ meet 
at a saddle point. Inspecting d'V, ii dp” reveals that p_ 
corresponds to the position of the local maximum of V,. In 
terms of £, the escape condition E = V,|,-, becomes 


a a 1/2 
e- [Zeue 3 +L) 2 = 9] . 2D 


2 
27 


Inverting this equation gives 
_ | 27€4 — 36€? + (9E? 
8(E€7 — 1) 


L£ 


3/2 1/2 
8) £] (23) 


The impact parameter b is defined as the perpendicular distance 
between the center of force and the incident velocity (Goldstein 
et al. 2001). It can be written as 

L L 


b=— = $ 24 
p ame (24) 


where P is the specific linear momentum. The maximum 
impact parameter for capture Dmax is given by 


[27€4 — 36€2 + (9€2 — 8)3/2€ + 8]! 


brax = -85 
24/2 (£2? — 1) 
The capture cross-section Ocap is given by 
2 m 27E4 — 36E? + (9E? — 8)3/7E + 8 
Ocap = Tbn = 5 5 : 
8 (er = 1) 
(26) 


Figures 1 and 2 are plots of bmax and the capture cross- 
section cap versus E, respectively. For ultra-relativistic 
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Figure 1. The maximum impact parameter for capture bmax vs. E for a neutral 
particle. 


Tcap 


Figure 2. The capture cross-section cap vs. E for a neutral particle. 


particles (€ >> 1), 


Dmax = 3v3 + V3 + o(=) (27) 
2 2E? E? 
The corresponding capture cross-section cap is therefore 
27r Or 1 
voy = 2 4 o( 5] 28 
P 4 28 E m 


For a slowly moving particle with speed v < 1, 


TARRA 29 
xi + —, 
5 (29) 
and thus 
2 2 
Dix = + OVE-1)= H OW), 30 
prt OWE=D => + 0H), GO 
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Figure 3. The value of £ at which p, and p_ meet (Lmin) vs. the 
electromagnetic coupling parameter a. 


and the capture cross-section becomes 


cap = z + Ov). (31) 
v 


4. Capture Cross-section of Charged Particles 


We will now follow the same procedure we used for the 
neutral particle. However, analytic expressions are not viable in 
this case and we will resort to numerical solutions, except in the 
ultra-relativistic particle case. The structure of the effective 
potential V, is generically similar to the neutral particle’s. The 
effect of œ is to raise (lower) the peak of V, for positive 
(negative) a. The effective potential V, is at an extremum 
when 


2af(p, — DE + o — LB — 2p.) — p. =0. G2) 


The extremum is a maximum when 
alL? — 2p.) + 3 — 4p.) 7] 
(p, DE + ph) px 


+ 2p, — 2£2<0. (33) 


To be consistent with the notation of the previous section, we 
let p, correspond to the minimum of V, and p_ correspond to 
the maximum. Here, Lmin (the value at which p_ and p, meet) 
depends on the value of a. The two parameters are related by 
the relation 


= +6001 (Cie — 3) — 8 La (Ca + 9) 


+ 3L hin (Chin — 3)? = 0. (34) 


Figure 3 is a plot of Lyin versus a. When a=1/2, Lmin 
approaches zero. This is because V, ceases to have a 
local minimum for a>1/2. Physically, this limit corres- 
ponds to the case when the Coulomb repulsion becomes too 
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Figure 4. The maximum impact parameter for capture bmax vs. € for a charged 
particle with a = 0 (black), a = —1 (blue), a = —2 (green), and a = —3 (red). 


Disc 


Figure 5. The maximum impact parameter for capture bmax vs. € for a charged 
particle with a = 0 (black), a = 0.1 (blue), a = 0.3 (green), and a = 0.5 (red). 


strong for stable orbits to exist as discussed in Al 
Zahrani (2021). 

Figure 4 shows how bmax depends on E for several negative 
values of the coupling parameter a. The effect of increasing |a| 
is to increase the values of bmax for all energies. This is 
expected because the Coulomb attraction makes it easier for a 
charged particle to get captured. In all cases, bmax is a 
monotonic function of E. In the ultra-relativistic limit, bmax 
approaches 3/3 /2, the limit in the neutral particle case, for 
any finite value of a, provided that a is not too large compared 
to €. 

Figure 5 shows how bmax depends on E for several values of 
a between 0 and 0.5. In this range, there is competition 
between the gravitational “attraction” and the Coulomb 
repulsion. The curves have richer structure. They fall quickly 
as E goes beyond 1 and reach a minimum. After that, the curves 
rise and reach 3./3 /2 asymptotically. 
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Figure 6. The max impact parameter for capture Dmax vs. E for a charged 
particle with with a = 0 (black), a = 1 (blue), a = 2 (green), and a = 3 (red). 


a 


Figure 7. The energy threshold for escape Emr vs. electromagnetic coupling 
parameter a. 


Figure 8. The capture cross-section Cecap vs. E for a charged particle with a = 0 
(black), a = —1 (blue), a = —2 (green), and a = —3 (red). 


Figure 6 shows how bmax depends on E for several positive 
values of a greater than 0.5. Generally, Dmax becomes smaller 
as qa increases. This is expected because the greater the 
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Figure 9. The capture cross-section Cecap vs. E for a charged particle with a = 0 
(black), œ = 0.1 (blue), a = 0.3 (green), and a = 0.5 (red). 
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Figure 10. The capture cross-section Gcap vs. E for a charged particle with 
a = 0 (black), a = 1 (blue), a = 2 (green), and a = 3 (red). 


Coulomb repulsion the more difficult it is for a charged particle 
to be captured. In fact, there is a threshold energy Emr below 
which capture cannot occur. It is given by 


Emr = at = (35) 
4a 
This equation is valid for a > 0.5 only. Figure 7 shows how 
Enr Varies with a. 

The capture cross-section Ccap corresponding to Figures 4, 5 
and 6 is shown in Figures 8, 10 and 9, respectively. In all cases, 
Ocap Versus E curves inherent the features of the Dmin versus E 
curves. 

For ultra-relativistic particles, we can write bmax as 


— 992 
bma = 343 _ 3a i 9 = 20% + (=) (36) 
2 E 6V3 E? E 
The corresponding capture cross-section is then 


2 
27x 9ra (4a +97 | of 1 ) (37) 
4 E 2E? 7 


Ocap = 
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These limiting results are in agreement with our numerical 
findings. 


5. Conclusion 


We have studied the capture cross-section of charged 
particles by a weakly charged Schwarzschild black hole. We 
have shown that a trace charge on the black hole can have 
prominent effects. 

When the Coulomb force between a charged particle and the 
black hole is attractive, it enlarges the capture cross-section 
significantly. This is expected since the Coulomb attraction 
enhances the capture of charged particles. However, when the 
Coulomb force between a charged particle and the black hole is 
repulsive, it shrinks the capture cross-section significantly. 
When the electromagnetic coupling strength is below a critical 
value, capture is possible for all values of the particle’s energy. 
When the electromagnetic coupling strength is above the 
critical value, there is a minimum value of the particle’s energy 
below which capture is impossible. This is because the 
Coulomb repulsion surpasses the gravitational attraction unless 
the particle’s radial momentum is large enough. 

Our results emphasize the assertion that charged black holes 
will favorably accrete charges of the opposite sign. However, it 
is still possible for the black hole charge to grow if the plunging 
charged particles are energetic enough to the limit that the 
capture cross-section becomes independent of the sign of the 
charges. Moreover, the fact that the electromagnetic coupling 
constant is three orders of magnitudes greater for electrons than 
protons suggests that it is relatively easier for a black hole to 
accumulate positive charge than negative charge. 

It will be astrophysically interesting to study the energies of 
charged particles near an astrophysical black hole to understand 
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better how the black hole’s charge evolves. The problem can be 
astrophysically more viable when other astrophysical black 
holes, such as rotating black holes, are studied (in progress). 
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